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SPACECRAFT SPIN AXIS REORIENTATION METHOD 

TECHNICAL FIELD 

[0001] The present invention is related generally 

to satellite reorientation systems. More 
particularly, the present invention is related to a 
system and method of reorienting the spin axis of a 
satellite . 

BACKGROUND OF THE INVENTION 

[0002] Spin axis reorientation maneuvers are 

utilized, over a spinning transfer orbit and during a 
mission, to position a satellite in a proper 
Perigee/Apogee burn attitude or a proper liquid 
apogee motor (LAM) burn attitude. A spin axis 
reorientation maneuver may also be utilized to change 
the sun polar angle of the satellite to provide 
favorable power and thermal conditions. The 
satellite spin axis can be reoriented via control of 
onboard thrusters . 

[0003] To perform a reorientation maneuver, a 

satellite controller compares a target attitude with 
an initial or estimated attitude to determine an 
angle error therebetween. The commanded and 

estimated attitudes are often represented by 
quaternions. The angular error signal can be 

represented by an error quaternion, or be referred to 
as the Euler angle by small angle approximation. The 
satellite then performs a closed-loop reorientation 
maneuver using the angular error to generate an 
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acceleration command signal is converted into torque 
generated by the thrusters and exerted on the 
satellite, which in effect moves the spin axis of the 
satellite toward a target orientation. Ability to 
precisely control the spin axis trajectory during the 
reorientation maneuver, or in effect, the movement of 
the satellite spin axis in inertial space, is highly 
desirable. Improved control of the spin axis 

trajectory can provide an increase in maneuver 
flexibility and accuracy, as well as simplify fault 
protection design and, as well as fuel efficiency. 
[0004] In particular, a ^minimum-angle" slew is 

often desired. Minimum-angle slew refers to the 
rotating of the spin axis about an axis that is 
normal to both a first unit vector, along the initial 
spin axis, and a second unit vector, along the target 
spin axis. In performing a minimum angle slew, the 
spin axis trajectory of the satellite follows an arc 
of a great circle, i.e., a generally non-curved, 
shortest distance path on along a perimeter of a 
sphere . 

[0005] Several reorientation control methods have 

been utilized to reorient a satellite with minimum- 
angle rotation. The control methods tend to generate 
a control signal, which is directed about an 
instantaneous eigenaxis or along the direction of an 
instantaneous angular position error vector, and 
include the use of reorientation and rate tracking, 
and angle and rate limits. The instantaneous angular 
position error vector is the vector representation of 
attitude error. However, these control methods are 
limited to being performed when the satellite is in a 
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non-spinning state and are incapable of controlling a 
spin axis trajectory during a reorientation maneuver. 
[0006] In a non-spinning state, the eigenaxis of 

the satellite is fixed in the satellite body 
coordinate frame and is stationary in inertial space. 
The eigenaxis is the Euler axis about which a single 
rotation can be performed to change the attitude of a 
body between orientations. In the spinning state, 
the instantaneous eigenaxis, in the body frame, is 
changing constantly. The constant change in 

instantaneous eigenaxis prevents the spin axis from 
following a minimum-angle reorientation trajectory 
using the stated traditional control methods of 
reorientation . 

[0007] Additionally, spin phase error ambiguity 

can negatively affect performance of a minimum-angle 
reorientation. An angular position error signal can 
be regarded as the composition of spin axis attitude 
error and spin phase error. The spin phase error 
refers to the angular error about the satellite spin 
axis. When the commanded and estimated attitudes are 
not synchronized in the spin phase, the direction of 
the position error vector is dependent on the 
magnitude of spin phase error. This is referred to 
as spin phase error ambiguity. The spin phase error 
ambiguity increases the difficulty in predicting the 
spin axis trajectory and maneuver time. 

[0008] Thus, there exists a need for an improved 

satellite reorientation system that allows for the 
accurate control of the spin axis trajectory and that 
allows for minimum-angle slew to be performed along a 
path of minimum distance end-to-end. 
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SUMMARY OF THE INVENTION 

[0009] The present invention provides spin axis 

reorientation systems and methods for a satellite. 
In multiple embodiments of the present invention, a 
reorientation controller for a satellite is provided 
that includes a slew rate command generator that 
generates a slew rate command signal ( a> r _ C md ) in 
response to an attitude error signal. The attitude 
error signal corresponds to the difference between an 
initial attitude and a target attitude. The slew 
rate command generator may introduce a phase lead (9 L ) 
into the slew rate command signal {a> r _cmd). The 
controller may perform a spin phase synchronization 
when the target attitude is unsynchronized in spin 
phase with the initial attitude. The satellite may 
perform a reorientation maneuver in response to the 
slew rate command signal (co r _cmd). 

[0010] The embodiments of the present invention 

provide several advantages . One such advantage that 
is provided by an embodiment of the present invention 
is the provision of introducing a phase lead into the 
generation of a slew rate command signal. The 
introduction of a phase lead can compensate for the 
slew rate response lag due to finite control 
bandwidth and control transport time delay. 

[0011] Another advantage that is provided by an 
embodiment of the present invention is the provision 
of performing spin phase synchronization when the 
target attitude and the initial or estimated attitude 
are not synchronized in spin phase. In so doing, the 
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stated embodiment avoids the uncertainty of spin axis 
trajectory caused by spin phase error ambiguity. 
[0012] Furthermore, another embodiment of the 

present invention provides spin axis trajectory 
shaping capability, which improves satellite 
reorientation control. This capability provides the 
advantage of increased flexibility of attitude 
maneuver . 

[0013] Moreover, the embodiments of the present 
invention in providing the above stated advantages 
minimize travel distance of a satellite during a 
minimum-angle slew of a spin axis reorientation 
maneuver. As such and in general, the stated 
advantages are more fuel efficient than techniques 
that involve moving the spin axis along an arbitrary 
path . 

[0014] The present invention itself, together with 

further objects and attendant advantages, will be 
best understood by reference to the following 
detailed description, taken in conjunction with the 
accompanying drawing . 

BRIEF DESCRIPTION OF THE DRAWINGS 

[0015] Figure 1 is a vector and block diagrammatic 

view of a satellite reorientation system for a 
satellite in accordance with an embodiment of the 
present invention; 

[0016] Figure 2 is a logic flow diagram 

illustrating a method of reorienting a satellite in 
accordance with multiple embodiments of the present 
invention; 
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[0017] Figure 3 is a signal flow and block 

diagrammatic representation of the reorientation 
system in accordance with an embodiment of the 
present invention; 

[0018] Figure 4 is a coordinate system 

illustrating the relationship between an 

instantaneous slew rate command signal and an 
estimated angular rate; 

[0019] Figure 5 is a two dimensional diagram 

representing the projection of a spin axis trajectory 
on to a plane normal to the satellite spin axis; 

[0020] Figure 6 is a coordinate system 

illustrating the introduction of phase lead in the 
generation of a slew rate command in accordance with 
an embodiment of the present invention; 

[0021] Figure 7A is a graph illustrating a 

simulated angular rate along an x-axis of a 
satellite, spinning at a first spin rate, for a spin 
axis reorientation maneuver in accordance with an 
embodiment of the present invention; 

[0022] Figure 7B is a graph illustrating a 

simulated angular rate along a y-axis of the 
satellite, simulated in Figure 7A, in accordance with 
an embodiment of the present invention; 

[0023] Figure 7C is a graph illustrating a 

simulated angular rate along a z-axis of the 
satellite, simulated in Figure 7A and 7B, in 
accordance with an embodiment of the present 
invention; 

[0024] Figure 7D is a spin axis trajectory plot of 

the satellite, simulated in Figure 7A-7C, in an earth 



6 



centered inertial coordinate frame in accordance with 
an embodiment of the present invention; 

[0025] Figure 8A is a graph illustrating a 

simulated angular rate along an x-axis of a 
satellite, spinning at a second spin rate, for 
another spin axis reorientation maneuver in 
accordance with another embodiment of the present 
invention; 

[0026] Figure 8B is a graph illustrating a 

simulated angular rate along a y-axis of the 
satellite, simulated in Figure 8A, in accordance with 
an embodiment of the present invention; 

[0027] Figure 8C is a graph illustrating a 

simulated angular rate along a z-axis of the 
satellite, simulated in Figures 8A and 8B, in 
accordance with another embodiment of the present 
invention; and 

[0028] Figure 8D is a spin axis trajectory plot of 

the satellite, simulated in Figures 8A-8C, in an 
earth centered inertial coordinate frame in 
accordance with another embodiment of the present 
invention . 





LIST OF SYMBOLS 


Aq 


Attitude error quaternion. 


A© 


Rotated position error vector. 


Ae lim 


Position error limit. 




Magnitude of the slew rate command. 


COr cmd 


Slew rate command signal . 


COscmd 


Spin rate command signal . 




Position feedback control gain. 
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Proportional rate feedback control gain 
matrix . 

Proportional rate feedback control gain. 
a c Acceleration command signal. 

cob Measured or estimated angular rate. 

Q)r Slew rate response. 

9 d Phase lag of the slew rate response co r . 

Oj Phase lead of a slew rate command signal. 

co s Magnitude of the spin rate offset command 

signal or spin speed of satellite. 
u 5 pin Unit vector along the spin axis. 

T delay Transport time delay of a satellite 

reorientation system . 
9 reor Reorientation maneuver angle. 

£ e Control position error signal. 

£ a Control rate error signal . 

DETAILED DESCRIPTION 
[0029] Spin axis excursion, which is defined as 



the deviation of the spin axis trajectory of a 
satellite from a non-curved path, has been determined 
to be due to finite control bandwidth, control 
transport time delay, and may also be due to spin 
phase error ambiguity. Finite control bandwidth and 
transport time delay cause the actual slew rate of 
the spin axis to lag in phase with respect to a slew 
rate command signal. This lag prevents the spin axis 
trajectory of a satellite to follow a minimum-angle 
slew for a spinning satellite during a reorientation 
maneuver. The embodiments of the present invention 
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account for this lag to provide an improved system 
and method of reorienting the spin axis of a 
satellite . 

[0030] While the present invention is described 
with respect to a system and method of reorienting 
the spin axis of a satellite, the present invention 
may be adapted to be used in various applications 
known in the art. The present invention may be 
applied in military and civilian applications. The 
present invention may be applied to aerospace 
systems, telecommunication systems, intelligent 
transportation systems, global positioning systems, 
and other systems known in the art . 

[0031] In the following description, various 

operating parameters and components are described for 
one constructed embodiment. These specific 

parameters and components are included as examples 
and are not meant to be limiting. 

[0032] Referring now to Figure 1, a vector and 

block diagrammatic view of a satellite reorientation 
system 10 for a spacecraft or satellite 12 in 
accordance with an embodiment of the present 
invention is shown. The satellite 12 is shown in an 
initial orientation 12A and a final orientation 12B 
of a reorientation maneuver. A first vector 16 
represents an initial orientation of the satellite 
spin axis 17 in inertial space. A second vector 18 
represents the target orientation of the spin axis 17 
in inertial space. Angle between the two vectors 16 
and 18 is referred to as the total reorientation 
maneuver angle 9 reor . The reorientation system 10 
maintains the spin rate a> 5 of the satellite 12 
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before, during, and after the reorientation maneuver. 
In one embodiment of the present invention, the 
center of mass movement of the satellite 12 along the 
orbit 21 is ignored, during the reorientation 
maneuver . 

[0033] The reorientation system 10 may include 

multiple sensors 14 for the estimation and 
determination of the attitude of the satellite 12. 
The commanded or target attitude may be determined 
and generated by a main controller 20, or received 
from another spacecraft (not shown) or ground 
stations 22 (only one is shown) via a transceiver 24 . 
The system 10 includes the controller 20 that 
performs a reorientation maneuver in response to the 
signals received from the sensors 14 and the signals 
regarding a desired attitude. The controller 20 in 
performing the reorientation maneuver generates a 
slew rate command signal co r _cmd . In response to the 
slew rate command signal co r ^cmd , the reorientation 
system 10 generates a torque signal, which is 
received by a control actuator 26 and used in 
repositioning the orientation of the satellite 12. 

[0034] The sensors 14 may include gyros 29, sun 

sensors 30, or other sensors known in the art for the 
estimation and determination of the attitude of the 
satellite 12. Any number of each of the sensors 14 
may be utilized. 

[0035] The controller 20 may be microprocessor 

based such as a computer having a central processing 
unit, memory (RAM and/or ROM) , and associated input 
and output buses. The controller 20 may be an 
application-specific integrated circuit or may be 
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formed of other logic devices known in the art . The 
controller 20 may be a portion of a central control 
unit or may be a stand-alone controller, as shown. 
[0036] The control actuator 26 is generally a 

torque-generating device and may be in various forms. 
The control actuator 26 may include multiple 
thrusters, actuators, motors, a combination of 
thrusters and reaction wheels, or may be in some 
other form known in the art . 

[0037] Referring now to Figure 2, a logic flow 

diagram illustrating a method of reorienting a 
spinning satellite in accordance with multiple 
embodiments of the present invention is shown. 

[0038] In step 100, the controller 20 compares an 

initial or estimated attitude of the satellite 12 and 
a desired or target attitude and generates an 
attitude error signal. The error signal may be 
represented by an error quaternion Aq. 

[0039] In step 102, the controller 20 determined 

whether spin phase synchronization is enabled. When 
the spin phase synchronization is enabled, step 104 
is performed, otherwise step 106 is performed. When 
the target attitude is synchronized in spin phase 
with the estimated attitude, the spin phase 
synchronization can be skipped or, in other words, 
step 102 may not be performed. 

[0040] In step 104, the spin phase synchronization 

is performed to remove the spin phase error from the 
error quaternion Aq, generated in step 100. Assume 
Uspin is a unit vector along the spin axis, the spin 
phase synchronization may be performed by the 
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quaternion multiplication Aq <— AqSq spin , where Sq spin is 
defined in Equations 1 and 2. 



The subscripts 1, 2, 3, denote the body-fixed control 
axes. The vector part of the quaternion Aq indicates 
the direction of the Euler axis. The scalar part of 
the quaternion Aq, or fourth component, is related to 
the rotation angle about the Euler axis. 
[0041] In step 106, an angular position error 

vector, represented by [A<f> , AG, Av|/] T , is generated 
based on the error quaternion. Where A<j> t AG, and Av|/ 
represent the error components about the three 
coordinate axes of the satellite body coordinate 
system. When the spin phase synchronization is 

enabled, the spin phase error is removed from the 
angular error signal in forming the position error 
vector [A0 , AG, A\y] T . 

[0042] Referring now also to Figure 3, a flow 

diagrammatic representation of the reorientation 
system 10 in accordance with an embodiment of the 
present invention is shown. 

[0043] In step 108, the controller 20 generates a 

slew rate command signal co r _cmd based on the position 
error vector [A<f> , AG, A\|/] T . 

[0044] In step 108A, a trajectory shaping logic 50 

of a slew rate command generator 52 rotates the 
position error vector [A<f> , AG, Av[/] T about the unit 

vector u S pi n by a constant angle to generate a rotated 




(1) 




e = Aqo)U spin ( i) + Aq(2)U spin (2) + Aqo)U spin (3) 



(2) 
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angular position error vector A0 . In one embodiment 
of the present invention, the trajectory shaping 
logic is a 3x3 rotating transformation or shaping 
matrix C sha ping/ which is described in further detail 
below. The slew rate generator 52, by multiplying 
the position error vector [A^ , A0, Ai[/] T by the 
shaping matrix C sha ping/ introduces a phase lead, ahead 
of the position error vector, into the slew rate 
command signal co r _cmd . 

[0045] Referring now to Figure 4, a coordinate 

system illustrating the relationship between the 
instantaneous slew rate command signal co r _cmd and a 

measured or estimated angular response a>b is shown. 
An arbitrary inertial coordinate system o-xyz is 
shown, whereby the spin axis of the satellite 12 is 

assumed to be initially along the vector oz 0 , and the 

target vector oz ( , represents the target orientation 
of the spin axis in inertial space. A first unit 

vector np\ is normal to a plane PI, which is spanned 
by the initial attitude z 0 and the target attitude z, . 
[0046] A minimum-angle reorientation maintains 

movement of the spin axis in the plane PI during the 
maneuver. However, when a slew rate command signal 

€O r _cmd is generated along a first unit vector n P \ , the 
slew rate response is actually along a second unit 

vector npi . The actual spin axis trajectory 60, as 
best seen in Figure 5, of the satellite 12 is tangent 
to a plane P2 . The fixed angle 0 d , between the 
planes PI and P2 , is referred to as the phase lag of 
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the slew rate response co r relative to the slew rate 
command signal co r _cmd - The phase lag occurs due to 
the finite control bandwidth and the transport time 
delay of the satellite 12 and can be determined using 
equation 3, where co s is the magnitude of the spin 
rate command or the spin speed of the satellite 12, 
and is the proportional rate control gain. 



0 d = arctan 



[0047] The first term arctan — is related to the 

control gain. T delay is the transport time delay, which 

is dependent on control system hardware and software 
implementations performed. The contribution of 

transport time delay to the phase lag 6 d is smaller 



than the first term arctan 



[0048] Without applying the trajectory shaping 

logic 50, which is described in further detail below, 
the spin axis trajectory undesirably follows a curved 
path, such as the trajectory 60. That is for any 
point Zi on the trajectory 60, the angle between plane 

P2 and line oz t is constant and is equal to the phase 

lag 6 d . 

[0049] Referring now to Figure 5, the coordinate 

axis Xi represents the remaining reorientation 
maneuver angle. At the initial attitude z 0 , Xi is 
equal to the total maneuver angle 6 reor . At the target 
attitude z t , Xi is equal to zero. x 2 represents the 
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spin axis excursion. The desired path, through a 
minimum-angle reorientation, is the straight line z Q z t 
from the initial attitude z 0 to the target attitude 

[0050] The parametric curve equations of spin axis 

trajectory in (x i# x 2 ) coordinates can be represented 
by equations 4 and 5 for attitude varying from 0° to 
90°. 

-_± 

*i =0 reor e^ cos</> (4) 
z£ 

The maximum spin axis excursion is equal to 
O reor e tm0d sin0 d , which is proportional to the phase lag 

[0051] The spin axis trajectory 60 can be changed 

by using a constant shaping matrix represented by 
equation 6. 

C shaping = Us P in Y + [/ 3 *3 ~ U *Pi» Y ] COS & L ~ [_ U spin x ] SHI 0 £ ( 6 ) 

7 3v3 is a 3-by-3 identity matrix, and u sp m x is the skew- 
symmetric matrix, which is represented and defined by 
equation 7. u sp i„ is a unit vector along the spin axis 
and represented in the spacecraft body frame. The 
spin axis may be any fixed axis in spacecraft body 
frame . 



Uspin* — 



0 Uspin^ Uspin™ 

Uspin^ 0 Uspin™ 

Uspin^ U spi „M 0 



(7) 
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[0052] The shaping matrix C shaping can be used to 

provide a phase lead 6 L on the slew rate command 

signal co r _cmd relative to position error vector [A<fi , 

AG, Av|;] T . By setting the phase lead 6 L equal to the 

phase lag 6 d , the controller 20 compensates for the 

phase lag in the slew rate response a> r and a minimum- 
angle slew can be achieved. 

[0053] Referring now to Figure 6, a coordinate 

system illustrating the introduction of the phase 
lead 0 L in the generation of the slew rate command 
signal co r _cmd in accordance with an embodiment of the 
present invention is shown. The lead angle 0 L is set 
equal to the phase lag 6 d . The planes PI and P2 are 
overlapping. Therefore, the resulting slew rate 
response co r is along the direction of the position 
error vector [&<j> , A9, Aiy] T . 

[0054] Although in the above -described method the 
phase lead 6 L is set equal to the phase lag B d , the 
spin axis trajectory 60 may also be shaped to a 
curved path on either side of the straight-line path 
by setting the phase lead 0 L to be less than or 
greater than the phase lag 6 d . 

[0055] Continuing on with the method steps of the 

embodiment of Figure 2, in step 108B, the slew rate 
command generator 52 compares the magnitude of the 

rotated angular position error vector A0 to a 
position error limit A0 lim , which is represented by 
switch 62 . 
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[0056] In step 108C, when the magnitude of the 

position error vector A0 is greater than the 
position error limit A0 lim a slew rate command signal 
co r _cmd is generated along the direction of the 

position error vector A0 with magnitude coi im . In 
step 108D, when the magnitude of the position error 

vector A0 is less than the position error limit A0 lim 

a slew rate command signal cOr^cmd is generated with a 
value of approximately zero . 

[0057] In step 109, when the magnitude of the 

position error vector A0 is greater than the 
position error limit A0 lim , the control position error 
signal e e is set approximately equal to zero. It is 

assumed that the position error limit A0 lim is small 
compared to the total maneuver angle 0 reor such that 
the control position error signal e e can be ignored 
during a large angle maneuver. 

[0058] In step 112, when the magnitude of the 

position error vector A0 is less than the value of 
the position error limit A0 lim the control position 
error signal s 0 is set approximately equal to the 

position error vector A0 . 

[0059] In step 114, the control rate error signal 

e a is generated. The rate error signal e m is 

defined as the difference between a desired angular 
rate and an estimated angular rate cob of the 
satellite 12. The slew rate command co r _cmd and the 
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spin rate command signal co & C md are summed followed by 
the generation of the rate error signal . The spin 

rate command signal a> s _ C md may be received from the 
ground stations 22. 

[0060] In step 116, the controller 20 generates an 

acceleration command signal a c through use of the 
position control logic and rate control logic. In an 
embodiment of the present invention, the position 
control logic is in the form of a position feedback 
control gain matrix K e and the rate control logic is 
in the form of a proportional rate gain matrix K a . 
The position error signal s e is multiplied by the 
feedback control gain matrix K e . The rate error 
signal e a is multiplied by the rate gain matrix K a . 
The proportional rate gain matrix K m is assumed to be 

a constant diagonal matrix with control gain k^ in 
the diagonal entries. The resulting products are 
summed to generate the acceleration command signal 

a c . 

[0061] In step 118, the controller 20 converts the 

acceleration command signal a c to a torque command 
signal. The control actuator receives and in 

response to the torque command signal adjusts the 
orientation of the satellite 12. 

[0062] The above-described steps are meant to be 

an illustrative example; the steps may be performed 
sequentially, simultaneously, synchronously or in a 
different order depending upon the application. 
[0063] Results of a first simulated satellite spin 

axis reorientation maneuver are illustrated in 
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Figures 7A-D in accordance with an embodiment of the 
present invention. For the stated embodiment, 

satellite spin rate was set at 4.5 degrees per second 
along the z-axis of the body coordinate system and 

the unit vector u sp m was set equal to [0 0 1] T . The 
shaping matrix C shap ing/ used for the stated simulation 
and the following simulation, is represented by 
Equation 8. 



C 

^shaping 



cos0 L sin0 L 0 
-sin0 L cos^ 0 
0 0 1 



(8) 



Satellite reorientation was initially controlled and 
balanced by reaction wheels. The reorientation 

maneuver arbitrarily started at 1000 seconds. The 
satellite angular rates cox, coy, and coz along x, y, 
and z axes are shown in Figures 7A-C. Figure 7D 
shows the trajectory 80 of the spin axis 17, or the 
trace of the unit vector u sp in , on a unit sphere in an 
earth centered inertial coordinate frame. The spin 
axis trajectory 80 closely follows a desired minimum- 
angle path 82 . 

[0064] Results of another simulated reorientation 

maneuver is illustrated in Figures 8A-D. The 
satellite spins at 30 degrees per second about a 
principal axis. The reorientation maneuver 

arbitrarily starts at 1000 seconds. The satellite 
angular rates cox, coy, and coz and the spin axis 
trajectory 90 are shown is Figure 8. Again the spin 
axis trajectory 90 closely follows a desired minimum- 
angle path 92 . 

[0065] The present invention provides a satellite 

reorientation system that incorporates the use of a 
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re -programmable shaping matrix to introduce phase 
lead into the reorientation of the spin axis of a 
satellite. The present invention also provides spin 
phase synchronization capability. As such, the 

present invention compensates for phase lag and 
removes spin phase error ambiguity, which in turn 
allows for accurate spin axis trajectory control of 
the satellite during a reorientation mission. 
[0066] While the invention has been described in 

connection with one or more embodiments, it is to be 
understood that the specific mechanisms and 
techniques which have been described are merely 
illustrative of the principles of the invention, 
numerous modifications may be made to the methods and 
apparatus described without departing from the spirit 
and scope of the invention as defined by the appended 
claims . 
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